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We compute the contribution to the Lagrangian from the leading order (2.5 post-Newtonian) 
radiation reaction and the quadrupolar gravitational waves emitted from a binary system using the 
effective field theory (EFT) approach of Goldberger and Rothstein. We use an initial value for- 
mulation of the underlying (quantum) framework to implement retarded boundary conditions and 
describe these real-time dissipative processes. We also demonstrate why the usual scattering formal- 
ism of quantum field theory inadequately accounts for these. The methods discussed here should be 
useful for deriving real-time quantities (including radiation reaction forces and gravitational wave 
emission) and hereditary terms in the post-Newtonian approximation (including memory, tail and 
other causal, history-dependent integrals) within the EFT approach. We also provide a consis- 
tent formulation of the radiation sector in the equivalent effective field theory approach of Kol and 
Smolkin. 

I. INTRODUCTION 

The anticipated detection of gravitational waves using ground-based laser interferometers [H, i, i, g has spurred 
significant advances towards a detailed theoretical description of binary compact objects. From a semi-analytical 
point of view, a major tool for understanding these systems comes from the post-Newtonian (PN) expansion, which 
is a perturbation theory based on the small relative velocity v of the binary's constituents. Chief among these 
accomplishments are the heroic calculations of the equations of motion for the orbital dynamics through 3.5PN order 
and the power emitted in gravitational waves, which are estimated to be sufficiently accurate for precisely estimating 
the parameters of the detected gravitational wave sources. See 0, Q for reviews of the standard treatments of the PN 
approximation and their main results. 

The recently introduced effective field theory (EFT) approach of Goldberger and Rothstein Q , called non-relativistic 
general relativity (NRGR), provides an economical way to apply the PN approximation to compact binaries. The 
NRGR approach is viable in any system wherein the wavelengths of gravitational perturbations are much larger than 
the sources or scattcrers. It is a valuable tool that utilizes the efficient diagrammatic and regularization techniques of 
perturbative quantum field theory to calculate quantities of interest within the post-Newtonian expansion. Only tree- 
level Feynman diagrams are relevant for (classical) astrophysical systems so that quantum corrections are completely 
ignorable. NRGR was further developed in Q and to incorporate spin effects and dissipation from gravitational 
wave absorption, respectively. It has also been recently used to derive the 2PN equations of motion for nonspinning 
compact objects in [1C| using Kol-Smolkin variables to streamline the number of contributing Feynman diagrams. 
NRGR has also been used in [l^] to derive the radiation reaction force on an electrically charged extended body 
interacting with its own electromagnetic radiation (see Section [II CI for discussion pertaining to the applicability and 
consistency of the particular method used in [H]). See [l^ for a good review of NRGR and also [M M, M, [la for 
applications beyond the two-body problem. In addition, NRGR has provided the first derivation of the 3PN spin-spin 
potential [H, [l^, US, HH and of a quadrupole-spin correction to the gravitational energy [23| . 

NRGR uses the language and formalism of perturbative quantum field theory via path integral methods. However, 
only the classical limit of the theory is relevant for describing astrophysical compact binary systems so that all quantum 
effects are blatantly ignored. Nevertheless, there are different path integral formalisms that can be used in NRGR to 
implement different boundary conditions on the gravitational perturbations. The boundary conditions are typically 
chosen to conform to the questions being asked. The "in-out" formalism, which is the framework typically discussed in 
many quantum field theory texts and was used in Q to introduce NRGR, is based on a scattering formalism that uses 
the Feynman Green's function to implement scattering boundary conditions on the gravitational perturbations at the 
asymptotic past and future. As such, the in-out construction is appropriate for calculating instantaneous conservative 
forces on the compact objects and for computing the instantaneous total power emitted in gravitational waves (at 
least, in a particular method), among other things. 
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One of the goals of this paper is to emphasize that the symmetric nature of the Feynman propagator, which is a 
consequence of the scattering boundary conditions employed in the in-out formalism, implies that it is unsuitable for 
self-consistently and systematically describing time- asymmetric processes related to dissipation and radiation reaction 
in compact binaries [53 • We stress that this is not a systematic flaw or shortcoming of NRGR but instead arises from 
not imposing retarded boundary conditions on the radiated gravitational perturbations. 

To address these and other issues in NRGR here we instead implement the "in-in" formalism to enforce retarded 
boundary conditions in a path integral framework. The in-in construction is an initial value formulation that evolves 
the system in real time from a given initial state and allows for the final state to be determined dynamically given only 
initial data. This is in contrast to the ab initio stipulation of the final state in the in-out construction for scattering 
processes. 

The in-in approach was first introduced by Schwinger psj as a way of computin g ex pectation values in quantum 
mechanics from a path integral formalism and was further developed by others in [2l.l25l. [26l . [27l [H, [2^, [13, HH, [H, Ull . 

Since its introduction, the in-in formalism has been extensively applied to problems where an initial value for- 
mulation is crucial for describing a system's dynamical evolution, typically involving nonequilibrium processes, from 
an initial state to an unknown final state. These include semi-classical gravity and stochastic gravity (see [34| and 
references therein), inflationary cosmology, quark-gluon plasmas, disoriented chiral condensates, thermal field theory, 
Bose-Einstein condensates and quantum Brownian motion, to name a few. See [ssf for corresponding references. The 
in-in formalism is also useful for addressing issues related to the quantum-to-classical transition (e.g., decoherence), 
macroscopic coherence, quantum kinetic theory, noise and fluctutations in open quantum systems, among other things 

In the extreme mass ratio inspiral scenario, the in-in formalism is crucial to guarantee the causal evolution of the 
binary in a curved background spacetime (e.g., Kerr) and has been successfully used in [3^ to rederive the first order 
self-force [H, [s^ acting on the small compact object. 

In this paper the in-in framework is used to derive, in the context of NRGR, the well-known compact object 
equations of motion with radiation reaction at 2.5PN order (first derived by Burke and Thorne [13, El, IH]) and 
the quadrupole gravitational waves emitted by the binary. The in-in formulation of NRGR should also be useful for 
deriving real-time quantities and the hereditary terms (e.g., memory and tail integrals) that appear in higher order 
expressions for the metric components (the gravitational waveform) and the radiated power. 

In Section [IT] we provide a brief description of the in-out formulation of the radiation sector in NRGR. In Section 
mil we provide a pedagogical presentation of the in-out framework and discuss its shortcomings for describing the 
real-time causal propagation of gravitational waves. We also review the in-in formalism. We then apply the in-in 
approach to derive the well-known results for the 2.5PN radiation reaction forces and for the emitted quadrupole 
gravitational radiation in Section IIVI In Appendix [Dl we apply the in-in framework to formulate the equivalent 
classical effective field theory (CIEFT) approach of Kol and Smolkin [i^ in a form suitable to self-consistently derive 
radiation reaction and other real-time quantities. This provides an alternative derivation of the in-in approach to 
NRGR. In this paper we focus on non-spinning compact objects and use the same conventions as 0- 



II. THE IN-OUT FORMULATION OF NRGR 



The central quantity in the NRGR paradigm is the effective action, Sgff- At the orbital scale r of the binary, 
the compact objects can effectively be treated as point particles interacting with nearly instantaneous potentials Hfj^^, 
and coupled to long wavelength (A ^ r), slowly varying, external radiation fields hf^,^. In the "in-out" path integral 
formulation the effective action is given by 



J J T^H^^ exp \^S[ri + /i + i?] + i ^ Spp[xK{t),ri + h + H^ 



(2.1) 



where S is the (gauge-fixed) Einstein-Hilbert action, Spp is the point particle action for each compact object in the 
binary and the index K labels each particle. 

Integrating out the potential gravitons i?^,^ from the theory at the orbital scale (in the Lorenz gauge and on the 
long wavelength background spacetime) schematically gives for the effective action 



/ 



Vhfj^u exp <.iS[rj + h] 




(2.2) 



where a curly line denotes a radiation graviton and the double solid line represents the compact binary. Wc take the 
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radiation graviton to be in the Lorenz gauge 



9"/ia/3 = ^dph\ (2.3) 

throughout the remainder. The diagrams in \2.2\ can be written perturbatively in powers of the relative velocity v of 
the binary by matching onto the effective theory at the orbital scale so that 



(2.4) 



Here L ~ mvr is the typical angular momentum of the binary. The first diagram in (j2.4p represents the leading order 
Newtonian contribution to the effective action, 



dt< 
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IX1-X2P 



(2.5) 



the second diagram provides the IPN corrections (from which the Einstcin-Infeld-Hoffman equations of motion [i^l 
are derived), and the third diagram gives the 2PN corrections, which have been recently computed using the NRGR 
approach in (lo| and agree with previous calculations using traditional methods (see 0)- 

The diagrams in (|2.4p do not couple to radiation gravitons, implying that there is no mechanism for the compact 
objects to dissipate energy (aside from possible tidal heating in neutron stars at higher orders) via radiation reaction. 
However, the remaining diagrams in (|2.2p describe the interactions of the compact binary with radiation gravitons 
and are responsible for driving the inspiral toward coalescence from the emission of gravitational waves. 

The diagrams containing a radiation graviton in (j2.2p represent a multipole expansion in terms of the size of the 
binary r and the wavelength of the radiation A. These lengths are related, since r/\ ^ v, which implies that these 
diagrams scale as powers of v. For example, the second diagram in (j2.2[) represents the full coupling of a single 
radiation graviton to the binary's center of mass and has a multipole expansion given schematically by 



= iSyl/2 + iSy3/2 £,1/2 + + ■ ' ■ (2.6) 



,1/2 



„3/2 



n5/2 



where the first three diagrams above are computed explicitly in Q, 

iS'^l/2^1/2 = 



iS, 



^3/2^1/2 



dt^mKhoQ{t,0) 

= / dt'S^ mfcx^T^hio(t,0) = in center of mass frame 

rripi I ^ 

i 



K 



2^1/2 



2771p; 



(2.7) 
(2.8) 
(2.9) 



and their scaling with (w, L) is given in the subscript. Here E is the leading order (Newtonian) total energy, L is the 
leading order angular momentum, Qij is the quadrupole moment of the binary 



rriK 



K=l 



1 



and i?oiOi is the linearized Ricmann tensor, given by 



^(^d^hij + didjhoQ - dodihjo - dodjhio^ + 0{h^). 



(2.10) 



(2.11) 



In obtaining these expressions, the radiation graviton field h^^{t,:ic) is expanded in multipoles around the center of 
mass, which is taken to be the origin of coordinates. Explicitly writing the velocity expansion into the effective action 
gives 



J Vh^i, exp|7S'pot[xi,X2,/i]|. 



(2.12) 



4 



„l/2 ,,5/2 „3/2 ^,3/2 



„9/2 „3/2 „7/2 



t,5/2 ,,5/2 „l/2 „4 ^1/2 

„7/2 „5/2 



uV2 „l/2 „l/2 „l/2 „3/2 „l/2 



„l/2 ,,3 „3/2 



FIG. 1: The diagrams that potentially contribute to radiation reaction at (a) 0.5PN, (b) 1.5PN and (c) 2.5PN orders, 
where 

Spot[^i,^2,h] = S[r] + h] + S^u"L^ +3^,21^1 + 5,u4j^i + h S'„i/2ii/2 + S'„3/2ii/2 + 51,5/2^1/2 H (2-13) 

is the (effective) action obtained by integrating out the potential gravitons. This expression for the in-out effective 
action is useful for deriving the conservative forces acting on the compact objects (tI. [lOj. for computing the instan- 
taneous power emitted in gravitational waves 0, Sll , and can even be used to recover the Schwarzschild metric of a 
single compact object [60j . However, to derive the radiation reaction forces in the equations of motion and to compute 
the multipole expansion of the emitted gravitational waves in NRGR requires imposing retarded boundary conditions 
by using an initial value (or "in-in" ) formulation to describe these quantities in real-time. To demonstrate that the 
in-out path integral formulation above is inadequate for these latter purposes, let us attempt to derive the leading 
order radiation reaction forces acting on the compact objects using (|2.12|) . 



A. Radiation reaction 



Integrating out the radiation gravitons involves calculating all tree-level and connected Feynman diagrams to the 
desired order in v and is equivalent to performing the path integral over h^i,. The resulting theory describes the 
motion of the compact objects (effectively as point particles) subject to radiation reaction due to the backreaction 
of gravitational wave emission. The diagrams that can contribute to radiation reaction through 2.5PN are given in 
Fig.©. 

In Appendix El we show that diagrams with a time-dependent (but otherwise generic) vertex coupled via a radiation 
graviton to a conserved quantity will not contribute to the equations of motion or the radiated power. The vertex 
associated with Syi/2j^i/2 is proportional to the total mass of the binary A/, which is conserved at leading order, 
implying that the single 0.5PN diagram in Fig.([T^) vanishes. Likewise, the v^^'^ vertex is zero in the center of mass 
frame; see (j2.8p . Therefore, both diagrams in Fig.lUJs) vanish. Similar reasoning implies that only the diagram 
containing two v^^^ vertices in Fig.(IT}:) does not automatically equal zero. Furthermore, the first two terms in (|2.9p 
are conserved at this order suggesting that these give a vanishing contribution in the effective action at 2.5PN. Indeed, 
only the quadrupole term in (|2.9p is relevant so that 

Fig.dlt) = iSlff"" = {2^) j '^^ J dt'Q,j{t){Ro,oj(t,0)Rokoi{t',0))Qki(t') (2.14) 

The leading factor of 1/2 is the symmetry factor of the diagram. The imaginary part of Fig.([TJ;) is originally computed 
in 0] and generalized to d spacetime dimensions in (45j . The real part of this diagram yields a contribution to the 
equations of motion upon variation with respect to the particle coordinates. 
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The two-point function of the product of Uncarizcd Ricmann tensors is 



{RoiOj{t, 0)Rokoi {t', 0)) 



20 



-ik"{t-t')_ 



,0\4 



fe",k 



(fcO) 



k2 



so that 



■ c2.5PN 



dt ( dt' Q,,{t)Q,,{t') [ e-^*^° — 



and the distributional identity in (jA3p imphes 
27116* AT 



■ qZ.SPN 



dt / dt'Q,,{t)Q,.j{t') / <^ PV 



fcO,k 



0\4p-ifc0(t-t') 



(fc0)2 - k2 + ie 



(2.15) 



(2.16) 



(fcO)4e 



(fc")2 - k2 



+ / (fc°)4e-"="(*-*')<5((/fc")2 - k^) 



(2.17) 



In the integral of the first term, t is constant so we may shift the variable to s = i' — t and expand Qij {t + s) in 
powers of s, giving 



•C2.5PN 



2TTiG 



N 
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1 d" 



dtQ^At)Y,-:^Quit) ds PV 

n=0 •' "'k 

(fc")4(5((fcO)2-k2)Qy(fcO)Oy(-fcO) 



'{ky 



kO (F)2-k2 



k,fcO 



(2.18) 



where Qij{k'^) is the Fourier transform of the binary's quadrupole moment; see (|A6|) for our conventions. In the first 
term, we recognize the master integral I{n, 4, 0) defined in (|Bip and evaluated in Appendix [Bj This implies that only 
the n = 5 term contributes, so that 



■ C2.5PN 
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d^Q,j{t) ttG 



N 



'knQ.,(iki)i^ 



{2n) 



(2.19) 



with 1(5,4,0) = 5!/(47r) from dig]). 

Our expression for the imaginary part of the effective action in (|2.19p agrees with 0, and yields the quadrupole 
formula for power loss upon inserting another factor of |k| to the integrand. The real part is proportional to the time 
integral of the Burke-Thorne "reaction potential" [i^ [4l|, 



(2.20) 



evaluated at the positions of the two particles. However, two integration by parts in the first term of (j2.19p reveals 
that 



RoSlfJ^ 



Gn 

To" 



dtQ^,{t)Q,At) = 



20 



2 

I] ) ! 



(2.21) 



which is a total time derivative and does not yield a contribution to the equations of motion. 

In the next section, we compute the emitted gravitational waves in the quadrupole approximation using the in-out 
effective action p.l2p and demonstrate that the gravitational radiation does not propagate causally from the source. 



B. Quadrupole radiation 

The first few diagrams contributing to the gravitational radiation emitted by the compact binary are shown in (j2.6p . 
In the physically relevant transverse-traceless (TT) gauge, the first two diagrams in (|2.6p vanish. This is expected 
since these correspond to monopole and dipole radiation, which are pure gauge effects and can be transformed away. 
This leaves the third diagram in (|2.6p . wherein the first two terms of p.9p do not contribute in the TT gauge. Applying 
the Feynman rules gives the leading order contribution 

1 = /i^-^(t,x) = A,j-fci(n) / dt' (/ife;(t,x)i?o,no«(t',0))Qmn(i') (2.22) 
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where n is the unit vector pointing in the direction of propagation and KijM is given expUcitly by 



Ajj,fci(n) = 6ik5ji - -^SijSki - njUiS^k - niUuSji + -nknAj + -n^rijSki + -niUjUkni, (2.23) 



which is pair- wise symmetric in ij, kl. See [4y| for a more detailed discussion regarding the relationship between the 
Lorenz and TT gauges. 

From the expression for the linearized Riemann tensor in (|2.1ip we see that the two-point function in (j2.22p is 

{hki{t, x)i?omOn(t', 0)) = ^(hki{t, x) \^d^,h^n{t\ 0) + d^djim{t' , 0) - dt'dmhnoit', 0) - dt'dnhmoit', 0)] ). (2.24) 
Only the first term contributes, so that 



^piHi (t, x) = ^A,,- fci(n) J dt' Dii^,^{t - i', x)Q,„„(t') (2.25) 



where we have integrated by parts twice to move the time derivatives onto the quadrupole moment. 



Taking the (physically relevant) real part of hj^ using 



Re*i?|;„„(i - <',x) = \(D]il^ + D^tn) = {{t - t'f - x^) P,.,,„„ (2.26) 

where Pkimn = {rjkm'nin + VknTiim - VkiVmn) and performing the time integral gives 

m-i'Rehjnt,^) = -^A,,,H(n) iQuit - |x|) + Quit + |x|)l . (2.27) 



We see that the gravitational waves do not causally propagate from the source because of the dependence on the 
quadrupole moment at the advanced time t + |x|. Furthermore, (|2.27p implies that there is no power emitted by the 
compact binary. 

To see this, recall that the power emitted in gravitational waves is the surface integral of the radial flux of gravita- 
tional waves t^^ [1^ 

P = J dAt'"- = - J dA {dohf^drhj^) (2.28) 

where angled brackets here denote averaging over many cycles of the gravitational wave and the surface integral is over 
a very large sphere centered on the binary's center of mass. Using (|2.27p and the relation drf{t± |x|) = ±(?o/(i ± |x|) 
one finds that 

iO'- (X (g,,(t- |x|)g,^.(i- |x|)) - (Q,^.(i + |x|)g,,(i + |x||)) (2.29) 

Both terms can be shown to have the same Fourier transform implying that the flux of purely retarded radiation is 
cancelled by the flux of purely advanced radiation leading to t^^ = and vanishing power loss. This explains the lack 
of radiation reaction in the equations of motion seen in the previous section. 

This raises an interesting observation about the in-out approach. Computing the power loss directly from the 
gravitational perturbations /i^;^ gives a vanishing result because the scattering boundary conditions imply that there 
is no net flux of gravitational waves leaving the compact binary in an irreversible and asymmetric manner. However, 
since the imaginary part of the effective action S^'^J^ is related to decay processes then a non-zero expression for the 
power loss can be computed in the manner discussed in Q. In particular, this quantity is related to the total number 
of gravitons emitted by the system during a large time T — > oo 

IW.,, - \ J ^-'^"^ (2-30) 

where dT is the differential rate for graviton emission. Therefore, the differential power is dP = uudT so that the 
integrated power spectrum is 

r d^r 

P = / dujdVluj—— (2.31) 
J diodil 

and was used to derive the famous quadrupole power loss formula in NRGR in 0, . 
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C. Radiation reaction, self-consistency and hereditary effects 



Before continuing, we mention some technical subtleties associated with the in-out calculation of the radiation 
reaction presented above. There we demonstrated that the initially divergent effective action in (|2.14p has a finite real 
part that is the integral of a time derivative, which implies that the radiation reaction force is zero when derived from 
a variational principle; sec (|2.2ip . However, the correct expression could be calculated by first varying the divergent 
effective action (|2.14[) and then regularizing the divergence in the resulting force [6l|. It is straightforward to show 
that doing so yields. 



x^^{t), (2.32) 



which is the correct radiation reaction force first given by Burke and Thorne [40|, |4l|, |42|. However, we remark that 
this procedure (of regularizing the variation of the action instead of the action itself) can lead to inconsistencies within 
the in-out formalism and to incorrect expressions for history-dependent or hereditary quantities. 

For example, the finite part of the real part of the effective action in (|2.I4[) gives the integral of a total time 
derivative (|2.2ip and yields a radiation reaction force equal to zero and not (|2.32p . Either one has a regular effective 
action that yields a vanishing radiation reaction or one has a divergent and ill-defined Lagrangian that provides a 
finite and correct expression for the radiation reaction force (at least, if it is local in time). Therefore, it seems that 
the in-out formalism does not give a finite effective action that simultaneously provides the correct expression for the 
radiation reaction force. 

In addition, the force (|2.32p is not conservative and represents a dissipative component to the compact object 
equations of motion that should correspond with the emission of gravitational waves satisfying retarded boundary 
conditions. However, the gravitational perturbations are still given by (|2.27[) and yield no power loss (when computed 
from the field directly as in (|2.28p instead of from the imaginary part of the effective action) despite the appearance 
of physical radiation reaction forces on the compact objects according to (|2.32p . 

Finally, consider the extreme mass ratio inspiral (see [47t for a comprehensive pedagogical review) . In this system, 
the radiation reaction (more appropriately called "self- force" ) is a force on the small compact object that depends 
on the entire past history of the small compact object's worldline and on the configuration of the gravitational 
perturbations it couples to in the past. As such, the self-force at proper time t is affected only by gravitational 
perturbations that have evolved from within the past null cone of the worldline. It is straightforward to show that 
using the in-out formalism to calculate the self-force by first varying the effective action and then regularizing the 
resulting self-force leads to an extra dependence from the future null cone due to the advanced propagator, thus 
violating causality. 

In general, the in-out formalism may not give correct expressions for quantities that ought to have a causal history- 
dependence. This likely includes the hereditary terms in the PN expansion of the near-zone metric, in the radiative 
field at future null infinity, and in the radiated power at higher orders (1.5PN relative order and beyond). If one 
wants to use NRGR to systematically and self-consistently compute the kinds of quantities discussed throughout this 
section then one should implement a framework different from the in-out prescription discussed here. 

In the next section, after presenting the underlying assumptions and construction of the in-out framework, we 
discuss the "in-in" approach, which provides a self-consistent formalism suitable for describing the causal evolution 
of the compact binary as it undergoes the real-time processes of dissipation and backreaction from the emission of 
gravitational radiation that causally propagates to future null infinity. 



III. COMPARING THE IN-OUT AND IN-IN FORMALISMS 



We strongly emphasize that the problems encountered in the previous section do not stem from a systematic flaw 
with NRGR. Rather, these issues follow from using a particular formalism for computing the effective action that 
docs not implement the (asymmetric) retarded boundary condition on radiation gravitons. This condition is necessary 
for self-consistently describing the real-time processes of dissipation (via radiation reaction) and backreaction (from 
gravitational wave emission) that occur during the dynamical evolution of the compact binary. In this section, 
we provide motivation for these statements and review the framework appropriate for addressing these issues. Our 
discussion will be somewhat pedagogical. The reader already familiar with the in-out and in-in formalisms can proceed 
ahead to Section HVl for the in-in representation of NRGR or to Section HV Al for the corresponding Feynman rules. 

We first present the in-out formulation of NRGR within the more standard notation of quantum field theory found 
in many common texts (e.g., [48l. [i^. [50| ) . This will allow for a more transparent presentation of the details regarding 
the in-out and in-in formalisms. 
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In the in-out path integral formulation (i.e., the one commonly encountered in textbooks), the generator of connected 
correlation functions is a functional of several auxiliary sources that couple to the radiation gravitons and to the 
coordinates of each particle's trajectory: 



j Y^VyiK J T>hf,^ exp\iSpot[:>Ci,y^2,h]+i^ J dt^K ■ + i J d'^x Jf^^h^^i , (3.1) 



where Spot is the (effective) action describing the two compact objects coupled to the radiation gravitons, which is 
computed by integrating out the potential gravitons from the theory 0]; see (|2.13p . Variation of W with respect to 
jfj^i^ gives the one-point function of the radiation graviton (i.e., the quadrupole radiation computed in Section [II B[) 



(3.2) 



while a similar variation with respect to ji<-(i) gives (xA'(i))in-out. A partial Legendre transform of W 



r[(xi), (X2), JH = W[h,h, Jn -J2 [ '^^j^ ■ (^if)in-out (3.3) 

K=l •' 

gives the functional whose variation of the real part yields the equations of motion for the orbital motion of the 
compact objects 



JRer 

= 



'5(xK(t))in-out 



(3.4) 

jl=j2 = ./""=0 



These equations describe the evolution of (x^(t))in_out: which are assumed to be the classical variables Xx(i) (since 
the masses of the compact objects are astronomically large implying that there are negligible quantum fluctuations 
ofiK(t)). 

Actually, (|3.ip describes a first principles description of the compact binary where the trajectories of the compact 
objects are quantum mechanical and the gravitons are quantum fields. As such, a mechanism for decohering the 
quantum trajectories Xi<-(i) is needed to ensure that the compact objects evolve within a (semi-)classical limit. 
However, this requires a more detailed discussion of the suppression of quantum interference effects, which is usually 
given in a density matrix formulation and is intrinsically an initial value problem. See [sH . [sil . [s^ . [s^ and references 
therein for further discussion on this and related points. 

Nevertheless, in the assumed classical limit, the path integral over the compact object trajectories is dropped and 
(xii-(t))in_out is taken as an external field from the point of view of the radiation gravitons. As we are about to show 
below, (xif(t)} in-out (as derived from W above) is not the expectation value of the position operator. Therefore, 
interpreting the (in-out) one-point function as a classical variable, which ought to be selected by physical mechanisms 
of decoherence, is inaccurate and can lead to certain problems when trying to describe the classical limit of a quantum 
theory. Indeed, classical variables are real and evolve causally in time. 



A. One-point functions in the in-out formalism 



We demonstrated earlier that using the in-out formalism does not give a contribution to the radiation reaction 
forces from an action principle and does not describe the causal propagation of gravitational waves emitted from the 
binary. In this section, we show why this occurs. 

Let us consider a much simpler theory to focus our attention on without having the distraction of cumbersome 
details that are otherwise irrelevant for this discussion. Consider a real, massive and linear scalar field in flat spacetime 
coupled linearly to a physical source (5(t, x) 

S[<P\ = \j d'^x [da,(j)d°'<i) + m^cf^ + 2Q(f\ . (3.5) 

The last term is analogous to the interaction terms, 5^1/2^1/2 which are linear in the radiation graviton field 

and depend on the trajectories of the compact objects. 
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In the interaction picture, a state |a) evolves from t = — oo according to \a;t) = Uj+Q{t, — oo)|q!) where J(t,x) is 
an auxihary source whose role will be obvious below. The time-evolution operator is 



Uj+Qit, t') = T cxp \ i dT J Sx [J(r, x) + g(T, x)] 07(r, x) I (3.6) 



where T is the time-ordering operator and (f)i is the field in the interaction picture, which is related to the Hcisenberg 
picture by 

07(i,x) = J7,7+Q(t, -oo)0H(i,x)J7j+Q(-cx),t). (3.7) 

The generating functional W[J] for connected correlation functions of the field operator is defined as [isl. 

^w[.j] ^ (o,out|J7j+Q(oo,-cx))|0,in). (3.8) 

Here |0, in) denotes the vacuum state in the distant past (the "in" vacuum) and |0,out) is the vacuum state in 
the distant future (the "out" vacuum), which does not necessarily equal |0,in), especially in a general background 
spacetime. The presence of these vacua in (|3.8p provides the name of this formalism, the so-called "in-out" approach. 
When the auxiliary source J and the interaction Q are adiabatically turned on and off (as is usually assumed in 
quantum field theory) we can set the in- and out-vacua equal to the state |0) since there is a preferred vacuum 
selected by the Poincare symmetry of the underlying background Minkowski spacetime. Therefore, we write (|3.8p as 



e 



^^[■^1 = (0|[/j+q(oo,-oo)|0). (3.9) 
The one-point function of the field is the functional derivative of W with respect to the source so that 

SW 



'^fl"(i,x))in-out 



(3.10) 

,7=0 



and is the analogue of the one-point function in (|3.2p . It is commonly implied that (|3.10p is an expectation value 
of the field because the in- and out-vacua are equal. We will now show that while the in-out one-point function (as 
derived from W) is the expectation value of some operator, it is not of the field. 
To see this, we use p.6p to show that 

^^^i±^|^^p^ =*f77+Q(cx),i)0/(t,x)J7j+Q(t,-^) (3.11) 

from which it follows that (|3.10p can be written as 

(0H(t,x))i„_out = {0\UQi+^,t)$i{t,^)UQit,-(x^)\0). (3.12) 
Writing the right side in terms of the Hcisenberg picture field, with ,7 = in p.7p . gives 

(0/7(t,x))i„_out = (O|J7Q((»,-«))0/f(t,x)|O) (3.13) 
- (g|0i/(t,x)|O) (3.14) 

where \Q) = J7q(— oo, oo)|0) is the state |0) evolved from the distant future to the distant past in the presence of 
the physical interaction Q. Notice that (|3.13p is not the expectation value of the field but of the operator Uq(I)h, 
instead. Therefore, one-point functions (and, more generally, ?7-point functions) in the in-out formalism are not true 
expectation values of the field. 

Instead, (|3.14p is a matrix element of the field operator for a non-vanishing interaction Q and is generally complex. 
For our real scalar field it is easy to show that p.l3p is not purely real, which also suggests that the in-out one-point 
function is not a classical field. (Recall that to make the graviton one-point function purely real, we simply took the 
real part of the generally complex expression in (|2.25p using (|2.26p .) 

Most importantly, and not unrelated, is the fact that the in-out one-point function docs not describe the retarded 
propagation of the field perturbations that are induced by Q. This can be seen by either computing the right side of 
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(|3.12p directly or using a path integral representation for e*^ and computing the variation as in (|3.10p . Both methods 
give 

(i,x))in-out = i j d^yDpix - y)Q{y) (3.15) 

= i y" d^y[Dret{x -y) + Dadvix -y)- iDnix - y)]Qiy) (3.16) 

after using (|C17p . Note the dependence on the advanced propagator and on the non-local Hadamard two-point 
function Dh, which is non-zero even when — y" is spacelike. While taking the real part of (|3.16p removes the 
latter, the advanced propagator remains and is responsible for the dependence at advanced times seen in (|2.27p . 

To summarize, the in-out one-point function, as derived from W in (|3.8p or p.9p . is actually a matrix element of the 
field operator between the (in-) vacuum and a state \Q) that depends on the details of the evolution of the interacting 
field (see the definition of |Q) below (|3.14p ). This matrix element is neither real nor causal and the final state of the 
system (|0,out} = |0} here) is specified ab initio instead of being determined dynamically from a given initial state, 
as would be appropriate when given only initial data. We will see below that a description of true expectation values 
is intimately related to an initial value formalism. 



B. In-in formalism for expectation values and causality 

To ensure that retarded boundary conditions are employed it is important to develop the appropriate generating 
functional. This is the so-called "in-in" generating functional, which was first introduced by Schwinger p3| and further 
developed by many others since (see Section U for references). Here we review the construction of the "in-in" effective 
action. 

Let the in-in one-point function equal the true expectation value of the field so that 

((^H(t,x))in-in = (O|0H(t,x)|O) (3.17) 

= {0\UQ{~oo,t)$i{t,x)UQ{t,-^)\0) (3.18) 

The first line implies that the one-point function on the left hand side is equivalent to the true expectation value of the 
field, and the last line shows that it is manifestly causal and real (Hcrmitian). The subscript "in-in" is included because 
the vacuum states that appear above are, more generally, the in- vacuum |0, in), which gives it the name "in-in." We 
seek to determine the functional whose derivative with respect to the auxiliary source yields these expressions. 
In analogy with (|3.9p . a first guess might be 



e 



■^[-'l = {0\Uj+q{^oo, oo)J7j+q(oo, -oo)|0). (3.19) 



However, p.l9p trivially equals one since the evolution operators are inverses of each other. Instead, we introduce 
an auxiliary source Ji that couples to the field during the forward evolution in time and a second source J2 for the 
evolution backward in time. 



e''^^-'''-''^ = (0|J7j,+q(-oo, ^)Uj,+q{^, -(X))|0). (3.20) 
Using (|3.1ip . the one-point functions are 

6W 



(0|^(i,x))i„_i; 



6W 



= {0\UQi^^,^)UQ{^,t)(j)i{t,K)UQ{t, -oo)\0) (3.21) 

Jl = J2=0 

= {Q\UQ{~(x,t)$i{t,x)UQ{t,^)UQ{oo,-oo)\Q) (3.22) 



<5J2(t,x) 

and equals the expectation value of the field in the vacuum state 

(0]j(t,x))i„_i„ = (,^2^(<,x))i„_i„ = (O|0^f(t,x)|O), (3.23) 

as required. 
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The functional in (|3.20p has a path integral representation given by 

gW[Ji,J2] ^ J x)^^2502 cxp i^iS[(l3i] - iS[(l32] + i J d'^xMi^i J d'^xJ2(l32^ (3.24) 

where = (/)2 at t = oo according to (|3.20p . This condition on the path integral contour provides an alternative 
name, the "closed-time-path" or CTP formalism. We will use "CTP" and "in-in" interchangeably throughout the 
remainder. Notice that the presence of two sources implies path integrals over two fields, one each for the forward 
and backward time evolution depicted in (|3.20|1 . 

For a linear scalar field, the path integrals can be evaluated exactly 



W[Ji, -^2] = 2 y d\ J d^y [Ja{x) + Qa{x)]G''^{x - ij)[JB{y) + Qsiy)], (3.25) 

which can also be computed directly from (|3.20[1 . Here A, B ~ {1, 2}, Qi = Q2 and 

= ( ) (3.26) 

where the quantities Gf,±,d are given in Appendix [Cl If Q itself depends on a dynamical quantity (as in the next 
section where the quadrupole moment depends on the particle coordinates) then Qi ^ Q2- 

To check that this generating functional gives the appropriate causal structure for the one-point function, we 
compute the variation of W with respect to one of the sources, using (|3.21[) or (|3.22[) . and after using again the 
identities in Appendix [Cl find 

(0ff(t,x))i„_i„= [ d^yGret{x-y)Q{y). (3.27) 



This demonstrates that the in-in one-point function manifestly satisfies retarded boundary conditions, which should 
be expected since p.20p depends on the initial state |0,in) — |0) only and not a final state ah initio. Indeed, the 
propagators G^^ always combine to ensure that quantities of physical interest, such as equations of motion and 
expectation values, are manifestly real and causal [23l. [26l. [3ll. [3^ . 



C. The Keldysh representation 



The two-point functions appearing in the matrix G"*^ are not all independent, implying that we can isolate the 
fundamental ones by a suitable change of basis [6^ . This is accomplished by writing the sources in terms of their 
average J+ and difference J_, 

(3.28) 
(3.29) 



J- = J1-J2 



J1 + J2), 



so that 



W[J+,J-]^'- J d^x J d^y[Ja{x)+Qa{x)]G-''{x-y)[My) + Qb{y)]. 
Here a, 5 = ± and the matrix of two-point functions is given in the so-called Keldysh representation p6j by 



^iGadv 
-iGret hGu 



(3.30) 



(3.31) 



where = and Gh is the symmetric (Hadamard) two-point function (see Appendix IC]) . The a, 6 indices are 

raised and lowered with the "metric" 



Cab = 



1 

1 



implying that J± ~ . The expectation value of the scalar field in the Keldysh basis is 

dW 



(0//(i,x))i„_i„ = (0H+(<,x))i 
and is easily shown to equal p.27p using (|3.30|) . 



lj±=o sj^ 



j±=o 



SW 



(3.32) 



(3.33) 



j±=o 
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IV. THE IN-IN FORMULATION OF NRGR 



The retarded propagation of gravitational waves and their backreaction on the compact binary in the form of radi- 
ation reaction can be self- consistently and systematically implemented using the initial value formulation of quantum 
field theory in terms of the in-in effective action introduced in the previous section. We show how this effective action 
gives rise to the (classical) equations of motion for the binary system and how the emitted gravitational waves (i.e., 
the gravitational waveform) can be computed using the in-in framework within NRGR. We begin by perturbatively 
expanding the effective action and writing down the corresponding Feynman rules that allows one to circumvent the 
explicit evaluation of the path integrals. 

Following (j3.24p a path integral representation for the in-in generating functional in NRGR is given by 



/ n 2?xkiPxk2Z^^i''Z^^2'' exp I iSpot[{yiKi}, hi] - iSpot[{yiK2}, h2] 

'' K=l I 

X! / ^^(j-'^i ■ ^Ki - jA'2 • Xif2) + i / d'^x{Ji^„h^{'' - J2^i.hi^'') } (4.1) 



where the conditions xa'i = 'X-K2 and /i^'' = h':^'^ are met at t = oo. (Gauge-fixing in the in-in formalism follows the 
well-known Faddeev-Popov [ssj procedure and is applied to both graviton fields, ft.^"' and h':^'^.) 

Perturbation theory in the in-in and in-out frameworks are similarly formulated. In particular, integrating out the 
radiation gravitons allows for (|4.ip to be written as 



^^W[{iKl},{iK2},Jl".J!i''] 



•' K=l [ K=l K=l'' 



+i / d'^xCi 



zo[jr,jn. (4.2) 



which is expressed as a certain functional derivative operator acting on a Gaussian functional of the sources and 
where the interaction Lagrangian is 

/CO 
d'xc^^t = {si:i[{^Ki},hi] - si:i[{^K2},h2]) ■ (4.3) 
n=l 

Here a superscript in parentheses denotes the number of radiation graviton fields contained in the interaction term. 
The quantity Zq is the free field generating functional for the radiation gravitons 

zo[jf ^ jn = / v^^v^" cxp _ ^^(2) ^ ^ j dix{ji^,h'^'' - j2^./ir) | (4.4) 

and is calculated by integrating the Gaussian along the CTP contour, which gives 

Zo[Jn = exp I - 1 Id^xl d'x' jf{x)Dfp.^,{x - xV^'W)^- (4.5) 
In the Keldysh representation, the matrix of free graviton two-point functions is 

Dl''^^,ix-x')=( , S^'>' (4.6) 

V ^-^aPfS 2^al3-iS / 

where a, 6 = ± and D^'^^g — 0. By construction, W yields true expectation values of both x^ and h^,y. 
Computing the (partial) Legendre transform of W gives the in-in effective action 

r[{(xK±)}, Jf] = W[{h<±}, - E / dffK ■ {^Ka) (4.7) 

K=l •' 

EE 5e//[{(iK±)},Jr]- (4.8) 



„5/2 j,5/2 



13 



FIG. 2: The non-zero diagram contributing to the leading order (2.5PN) radiation reaction in NRGR using the in-in formahsm. 
Note the CTP indices a, b that relate to the forward and backward branches of time in the closed-time-path contour of the 
CTP path integral. 

The equations of motion for the true expectation values of the particle coordinates are then found by varying the 
effective action, 

(4.9) 

One can show that the large masses of the compact objects and the influence from graviton quantum fluctuations 
sufficiently decohere the coordinate operators x^f and ensure that the above equations of motion for the expectation 
value appropriately describe the classical limit of the compact binary. See [5l|, [13, IH, HI] and references therein for 
elaboration of this point. The gravitational waves radiated by the binary can be computed from the graviton in-in 
one-point function (see (|3.33ll ) 

(4.10) 

X/c-=0, XK+=:x.K,iK± = J±"=0 

and forms the vacuum expectation value of the radiation graviton field. 

One can generally show that (|4.9p yields equations of motion for the trajectories of the compact objects that are 
manifestly real and causal. Furthermore, (|4.10p is a solution to the field equations for the radiation modes and is 
guaranteed to satisfy retarded boundary conditions. See [s^, HI] for proofs of these statements. See also [3^ for 
a more thorough discussion of the in-in/CTP formalism and its use for describing dissipative and, more generally, 
nonequilibrium processes in quantum field theory and quantum mechanics. 

A. NRGR Feynman rules in the in-in formalism 

The diagrammatic structure of in-in perturbation theory is nearly identical to the in-out approach. However, when 
drawing Feynman diagrams in the in-in formalism it is necessary to include CTP labels a, 6, c, . . . at each vertex to 
keep track of the forward and backward branches of time in the CTP path integral. For example, the leading order 
contribution to radiation reaction in the in-in approach is given in Fig. ([2]). 

The Feynman rules for calculating the in-in effective action (e.g., for deriving radiation reaction forces) are similar 
to those given in Q: 

• Include a factor of the radiation graviton two-point function D^^^^,g,{x,x') connecting particle-field vertices 
labeled by CTP indices a and b at spacetime points x and x' , 

• A vertex coupled to n radiation gravitons is labelled with n CTP indices, one for each graviton, 

• Sum over all CTP indices. 

When computing the one-point function of the radiation graviton there is an additional rule for those graviton lines 
that connect to a distant field point. In the Keldysh representation: 

• Include a factor of D~p^igi{x, x') for each radiation graviton that connects a vertex at event x' with CTP index 
a to the field point x. 

In the following sections, we use these rules to derive the expressions for the leading order radiation reaction and 
gravitational wave emission. Initially, we use the matrix of propagators given in (|4.6p to demonstrate these rules. 
However, one can cquivalently use (|4.19p instead for all tree-level diagrams since (|4.19p contains only the relevant 
propagators for constructing the effective action and graviton 1-point functions in the classical limit as we show below. 



5T 



1% ( W - 
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B. Radiation reaction 



The diagram for the leading order radiation reaction is given in Fig. ([2]). From the Feynman rules above we write 
down the term in the 2.5PN effective action that corresponds to Fig. ((21), 



2m 



pi 



dt / dt'Qjf (0(i?o% Wi?U(i'))Q''(i') 



(4.11) 



where we are using the summation convention for the CTP indices a,b = ±. The two-point function of the product of 
linearized Riemann tensors is as in (|2.15p but replacing the Feynman propagator there by the in-in two-point functions 
KB so that 



where D-^,^ 



PijkiD'''' defines D"'' and 



SikSji + SikSji - -6.ij6ki 



1 



in four spacetime dimensions. Performing the CTP summations gives 
1 



(4.12) 



(4.13) 



dt / dt' 



-2^Qm^^Dret{t - t', mtit') + \Qm^4oDH{t - t', Q)Q^l{t') 



' dt^ dt'^ 



where the differenced and averaged quadrupole moments are 

Q^lit) ^ Q\^t)^Q'^{t) 



nit) ^ l{Qnt)+Qi'{t) 



(4.14) 

(4.15) 
(4.16) 



and Q^i2 a-re the quadrupole moments for the worldline histories {x^^}. In terms of the xa'± variables one can show 
that 



'it) 



(4.17) 
(4.18) 



According to (|4.9p we can ignore the last term in (|4.14p and the 0{x^) terms from Q!^ since these do not contribute 
to the equations of motion. In fact, the 0(x?_) terms are related to the effects that the quantum fluctuations of the 
radiation graviton field have on the trajectories. This is discussed in more detail for the case of a particle moving in 
a curved background spacetime in [s^ . Issj . We will ignore these higher order terms throughout the remainder since 
they do not contribute to the equations of motion or the gravitational wave emission in the classical limit. In fact, in 
the Feynman rules we can (and from here on we do) simply use the following matrix of propagators 



^aP^six x) — ( -ryret 



■ piadv 



(4.19) 



instead of 

The 2.5PN effective action can then be written as 



zS!fP[{^K±}] = dt J dt' ^-^Dret{t-t\0)Y,mKx),_it)x^^4t)Q%it') 



where we ignore the 0(x^) terms in Q^f. Using 

Dret{t-t',0) 



2-K 



e{t - t')5{{t ~ t') 



l\2\ 



(4.20) 



(4.21) 
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and dimensional rcgularization one can show that 



from which the effective action follows 



•' K=l 



and provides the following term to the equations of motion 



rc'2.5PN 



5 dt^ ""^^^^ 



(4.22) 



(4.23) 



(4.24) 



with ; 



327rG'Ar. We observe that this is precisely the radiation reaction force on the K body obtained by Burke 



and Thorne [40|, l4l|, 143, \4t 



C. Quadrupole radiation 

The diagram for the leading order (quadrupole) radiation is given in (|2.6p with the addition of the CTP labels 
and "a" at the free end of the radiation graviton and at the vertex, respectively, so that 



(4.25) 

X/f-=0, X/f+=XK 



in the TT gauge. The factor of 1/2 in the square brackets comes from the linearized Riemann tensor. Only the a = + 
term contributes in the CTP summation so that 

m-^hj^it^ X) = -^K,.ki I dt' DltJt - t', x)Q""(t'), (4.26) 



with 



Dlt^{t-t',^) = -eit^t')^-^^L^^-^P,,„,„. (4.27) 

It follows that the quadrupole radiation is causally propagated from the source, 

m;/^^(t,x) = -^A,,,«QfeKi - |x|). (4.28) 

Notice that (|4.28p is purely real and causal, which is guaranteed by the in-in formalism, and that the causal and 
time-asymmetric emission of gravitational waves is consistent with the presence of radiation reaction forces in (|4.24p . 

From (|4.28p one can compute the leading order contribution in a multipole expansion to the power radiated per 
unit solid angle 



\^(k[^(t,^h^{t,^) (4.29) 



dP 

dn 

from which the total power is given by 

P-^(Q^,it~\^\)Q.,{t-\^\)) + --- (4.30) 

where here the angled brackets denote a temporal average over several periods of the gravitational wave and • • • 
represents higher order contributions in the multipole expansion. Notice that the power depends on the quadrupole 
moment of the binary at retarded times. 

Interestingly, the multipole expansion of the power loss and change in time of the linear and angular momenta 
can also be computed directly in the in-in approach using Feynman diagrams. Indeed, the expectation value of the 
graviton stress tensor is related to the coincidence limit of the Hadamard two-point function, D^p^g{x, x') [s^ . and 
is in turn related to the imaginary part of the effective action in (|4.14p . 
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FIG. 3: The coupling of a radiation graviton to a time-independent vertex C"'^ and a general, time-dependent vertex V°''^(t), 
represented here as a gray circle. 

V. CONCLUSIONS 

In this paper we computed the leading order radiation reaction forces on the compact objects within NRGR using an 
initial value formulation of the underlying path integral framework to implement retarded boundary conditions on the 
radiated gravitational perturbations. We also calculated the quadrupole radiation emitted causally by the compact 
binary and the leading order contribution to the radiated power. We showed that using the in-out framework is not 
suitable for describing the retarded propagation of true expectation values of radiation gravitons and their backrcaction 
on the compact objects in the form of radiation reaction. Using the in-in formalism guarantees real and causal n-point 
functions and equations of motion for expectation values. 

Whereas the in-out formalism is usefully applicable for the practical purposes of describing the conservative forces 
appearing in the equations of motion for the compact objects and the instantaneous power emitted in gravitational 
waves, the in-in formalism is necessary for self-consistently deriving the radiation reaction forces on the compact objects 
and the multipolc expansion of gravitational waves propagating causally in the far-zone, among other things. The 
self-consistent nature of the in-in approach to NRGR should also be useful for computing the hereditary contributions 
to the power loss and radiation field. 

In Appendix ID] we discussed the equivalence of NRGR and the classical effective field theory (CIEFT) approach of 
Kol and Smolkin [4^. We also showed that using only retarded propagators does not lead to a consistent theory in 
the original formulation of [43j . We provided a new formulation of CIEFT, which is simply an alternative derivation 
of the in-in approach to NRGR, that is appropriate for the radiation sector and consistently implements retarded 
boundary conditions within a strictly classical framework. This is achieved by doubling the degrees of freedom in a 
manner that is equivalent to the in-in formulation of NRGR. The appearance of advanced propagators at intermediate 
steps is necessary to ensure the self-consistency and causality of the derived equations of motion and gravitational 
waves. 
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APPENDIX A: COUPLING TO CONSERVED QUANTITIES 

In this Appendix we show that the coupling of a generic vertex V'^'^{t) to a conserved quantity C"^ gives a vanishing 
contribution to the equations of motion. Applying the (in-out) Feynman rules to the diagram in Fig.Q gives 

/OC />OC 
dt dt'iC'''''D^p^sit-t',0)iV\t'), (Al) 
-oo J —oo 

where C"-^ is time-independent to the given order in velocity. The vertex V^^lt') could be built out of other radiation 
gravitons but its exact structure is irrelevant for this general discussion. The momentum space representation for the 
Feynman propagator gives 



Using the distributional identity [4£ 
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where Py denotes the principal value of the integral and f{k'^) is an arbitrary function, yields an expression more 
amenable for isolating the real and imaginary parts of S 



k 



In the first term, define s = t' — t and expand ¥^^^[1 + s) in powers of s so that 



(A5) 

where the Fourier transform of V^^{t') is 

Vt'''{k°)= dt' e'''°*'V'"'{t'). (A6) 



From the discussion in Appendix [B] the only contribution from the master integral to the first term in (|A5|) occurs 
when n ~ 1. Also, the second term vanishes entirely since 

dt f f e-'''"'S{{k°f -k^)V^\k°) ^V''\0) f Sik^) ^0. (A7) 



k JfcO 

Therefore, the contribution to the effective action from a diagram with the structure given in Fig. ([3]) is 

dt—^, (A8) 

which is a total time derivative and does not contribute to the equations of motion, as claimed. It is a straightforward 
matter to extend the calculation here to the in-in case. 

APPENDIX B: MASTER INTEGRAL 

The master integral we use in the main body of the text is defined as 

/CO r p ikP s 

and is non-zero when q is an even integer. Writing s" as n derivatives of e''^ ^ with respect to then integrating 
over s and gives 



Iin,p,q) = i 



ld(fc") 



(B2) 

k'>=0 



Under the k integral the product of momentum vector components can be written as 



|k|9 

lk|9 



- 7^7\^{^n^.■■■S^,^^^,+■■■) (B3) 
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where the quantity in parentheses represents all possible pairings of the q indices [46[. In spherical coordinates, the k 
integral becomes 

/oo jd-lL 1 r poo 

and integrating over the d — 2 dimensional sphere gives 



I{n,p, q) = i 



(q + l)!!2rf-27r^r(i^) \ d{k") 



poo ,d-2+q 1 ~| 



(B6) 

A;0=0 



Using dimensional regularization to evaluate the remaining integral and computing the fc*^ derivatives on the resulting 
expression yields 



2d-i7r— r(£bl) (g + l)!! (p + g + d-3-n)! 



Notice that there is a non-zero contribution only when p + q + d — 3 = n implying that the master integral in d 
spacetime dimensions is 

Iip + q + d-3,p,q)= ^"'''^'^ec(f^) , (gg) 



and in = 4 dimensions becomes 



Iip + q+l,p,q)^-^--^S,^...,^. (B9) 
For n^p + q + d — 3, the master integral vanishes. 

APPENDIX C: PROPAGATORS AND TWO-POINT FUNCTIONS 

In this Appendix we collect some definitions, identities and relations for the quantum two-point functions that are 
relevant for this work. 

The positive and negative frequency Wightman functions are defined as 

^tp-y'S' (2;, x') = {hapix)hy>S' (x')) (CI) 

^a/iys' (x, x') = (hys' {x')Ki3{x)) , (C2) 
respectively. The Feynman, Dyson, Hadamard and commutator two-point functions are, respectively, 

Dai3t'S'{x,x') = {Thap{x)hys'{x')) (C3) 

Df,p^,s,{x,x') = {T*h^p{x)hys'{x')) (C4) 

D%^,Ax,x') = {{Kp{x);h,,s'{x')}) (C5) 

= {[Kp{x):h^,s'{x')]) (C6) 

where T is the time-ordering operator and T* is the anti-time-ordering operator. The field commutator is independent 
of the particular state used to evaluate it. Given the Wightman functions in (jCip and ()C2p we write the above two- 
point functions in the form (ignoring the tensor indices from here on) 

Df{x,x') = e{t - t')D+{x,x') + e{t' - t)D^{x,x') (C7) 

DDix,x') = 9{t' - t)D+{x, x') + e{t ~ t')D_{x,x') (C8) 

Dh{x,x') = D+{x,x') + D^{x,x') (C9) 

Dc{x,x') = D+{x,x') - D^{x,x'). (CIO) 
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From these we define the retarded and advanced propagators by 



- iDretix, X') = e{t - t')Dc{x, x') (Cll) 

+tDadvix, x') = e{t' - t)Dc{x, x'). (C12) 

These propagators also satisfy tlie following useful identities 

-iDret{x,x') = Df{x,x') - D^{x,x') (C13) 

= D+{x,x')^Dd{x,x') (C14) 

iDadv{x,x') = Dd{x,x) - D^{x,x) (C15) 

= D+{x,x')-Df{x,x') (C16) 

from which the Feynman propagator can be written in terms of its real and imaginary parts as 

Df{x,x') = ^'-[Dret{x,x')+Dadv{x,x')\-]^DH{X,X'). (C17) 

The Feynman. Dyson and Wightman functions are not all independent since 

Dh{x,x') = Df{x,x') + Dd{x,x') (C18) 

= D+{x,x') + D_{x,x'). (C19) 

Under the interchange of x and x' the Feynman, Dyson and Hadamard two-point functions are symmetric, the 
commutator is anti-symmetric and 

D+{x,x') = D-{x',x) (C20) 

Dret{x,x') = Dadvix',x). (C21) 



APPENDIX D: GRAVITATIONAL WAVES IN CLASSICAL EFFECTIVE FIELD THEORY 

The classical effective field theory (CIEFT) approach, introduced in [i^ and developed further in [Til. IZsl. IstI. Is^ . 
is an attempt to simplify NRGR [4J| by using only classical methods and techniques. More specifically, CIEFT uses 
Newton's constant Gat in the power counting instead of the Planck mass m^i^ = 32TrGN, vertices and propagators are 
defined without factors of i, integrating out gravitons from the theory does not use non-rigorous path integral methods, 
and radiation graviton lines in Feynman diagrams represent retarded propagators since one is working within a classical 
field theory with retarded boundary conditions [G^I • The first two differences are superficial since these are related to 
subjective conventions [6^ . The third attempts to eliminate the use of path integrals and to integrate out gravitons 
by solving for the short wavelength perturbations on the long wavelength classical background, then evaluating the 
resulting action [i^. Of course, this is equivalent to integrating out gravitons in NRGR by using the saddle point 
approximation to evaluate the path integral, viz., by drawing all tree- level Feynman diagrams. Calculating potentials 
and power loss in both NRGR and CIEFT boils down to the same set of Feynman rules so using path integrals (which 
are not directly evaluated anyway) or not does not seem like a significant simplification or modification of NRGR. 
Finally, the fourth embodies an attempt to apply retarded boundary conditions to the radiation sector of NRGR 
using only classical methods and quantitites. At the orbital scale, perturbations propagate (nearly) instantaneously 
so there is no difference in using the Feynman or retarded propagators for calculating potentials. Hence, CIEFT is 
equivalent to NRGR at this scale since both approaches describe the same system in the same classical limit using the 
same techniques (Feynman diagrams and dimensional regularization) with the same Feynman rules and propagators. 
To avoid confusing the reader, we emphasize that in computations to date, CIEFT and NRGR are equivalent and we 
are reluctant to distinguish CIEFT as an approach that is markedly different from or simpler than NRGR. 

However, it is worth pointing out that implementing retarded boundary conditions in a purely classical framework is 
more involved than stipulating that radiation gravitons represent retarded propagators in Feynman diagrams. Indeed, 
this Feynman rule is not correct as it leads to an inconsistency. 

To see this, we compute the quadrupole radiation. From the third diagram in (|2.6p this is found in the TT gauge 
in CIEFT to be 

- 1 r 2G 

m^iX^ = -^^^vM j dt' Dlturit - t', x)g,„„(t') = -^A,^- ,;Q,;(t _ |x|), (Dl) 
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which is the correct expression. However, the quadrupole radiation can be computed directly from the effective action 
by coupHng the radiation graviton to an auxihary source J^^ (x" ) and taking the variation of Sef / so that 



(D2) 



The effective action is calculated by solving the wave equation subject to retarded boundary conditions. Reinserting 
this solution back into the action, which is what is meant by "integrating out" [i^ , it is easy to show that 

Seff -Ifd^xJ d^x' [S^^'^ix) + J^%x)]Dl^:^,{x - x')[S-P{x') + r^{x')] (D3) 

where S^y represents the quadrupole source for the gravitational perturbations and is defined by 

V(t,x);i^''(t,x) = ^!iWi?o«J-(t,0),53(x). (D4) 

(S^'^ should be considered as an operator since the linearized Riemann tensor contains two derivatives acting on 
the graviton field.) We assume that the radiation gravitons couple to the auxiliary source bilinearly, / d^x J'"'/ip,y. 
Relabeling the integration variables t ^ t' allows us to write the retarded propagator in (|D3[) as the average of the 
retarded and advanced propagators. Taking the functional derivative of (jD3p . as in (|D2p . gives in the TT gauge 



:iW(i^'') = ^^^vM I dt' 



-"'-pi J 

which is not equal to (jDl[) . 

The Feynman rule for radiation gravitons in CIEFT is therefore seen to be 



Qrnnit'), (D5) 



(11^ -0 0000000 - = ^ {D;%ix ~y) + DfZ^ix ~ y)) (D6) 

and is not given by the retarded propagator alone. Indeed, this rule is equivalently given by ^^[iDj^^^p] and is 
therefore closely related to the in-out formalism of NRGR. As such, it is easy to show that there is no corresponding 
radiation reaction term in the equations of motion. The problem lies with the fact that the EFT paradigm works 
at the level of the action and not the equations of motion. Implementing retarded boundary conditions at the level 
of the Lagrangian requires something different from demanding that all gravitons be retarded solutions to the wave 
equation. 

Recall that imposing retarded boundary conditions on the radiated perturbations in NRGR effectively amounts 
to doubling the degrees of freedom such that each set of variables combines in the appropriate way to ensure that 
the binary's equations of motion and radiated gravitational waves are manifestly real and causal. Motivated by the 
in-in formalism, we generalize the radiation sector of CIEFT so that radiation reaction and gravitational waves are 
appropriately derived within a purely classical framework. Indeed, this may be viewed as an alternative derivation of 
the in-in approach to NRGR. 

To make the presentation below clearer we ignore non-linear interactions in the radiation sector. However, their 
inclusion does not introduce any conceptual obstacles to the framework. 

Consider the action for the compact objects and gravitational perturbations that results from doubling the degrees 
of freedom in a manner analogous to the in-in construction 

S[{^Kl}, {XK2}, /if^] = S[{^Kl}, K] - S[{^K2}, K"] m 

where, in the Lorenz gauge, 

S[{^ka}, h^] = ^ / (df^hAapd'^hf - Id^hAd^'hA^ + J d^x{jAc.p + SAc.p)hf + 0{h^), (D8) 

h is the trace of the radiation graviton field and A = 1,2. Here Jfj are auxiliary sources and denotes the 
quadrupole source as before. The first term in (jD7p is the usual action from classical field theory while the second 
term can be regarded as the action for the system (with 1-^2) evolving backward in time. We remark that (jD7p is 
the same action appearing in the CTP path integral of (|4.ip with ]ki = jK2 — 0. Also, (jD7p is not an action typically 
considered in classical field theory but seems necessary to achieve our purposes here. 
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Working with the h^!^ variables gives 

S[{^K±}, m = ^ / {d.hl^d^h:^' - Id.h'^d'^ha^ + J d^xiJ^f, + SlpYK^ + 0{hl) (D9) 

where the metric Cab is used as before to contract the CTP indices. Integrating out the gravitational perturbations 
involves first solving the wave equations 

P^'^pUhf = J^" + 5f + 0{hl) (DIO) 



subject to the appropriate boundary conditions and inserting the solutions back into (jD9p . In other words, we 
extrcmize the action with respect to variations in the radiation graviton fields and then compute the extremal "value" 
of S to obtain the effective action. 

As with the in-in formalism, we will set xa'i = ^K2 = and Jf^"^ = at the end of all functional variations, 
which implies that, momentarily dropping the spacetime indices, h and /i_ — ^- 0. Notice that the equation 

is sourced by terms that remain non-zero when x/^i = Since /i+ corresponds to radiated gravitational waves in 

this limit then these must satisfy retarded boundary conditions so that 



h+,A^)= J d''yD;i'^^{x~y)[Sfiy) + jfiy)]. (Dll) 
Solving the /i_ equation in (jD10[) under momentarily unspecified boundary conditions gives 

h_^,{x) = f d^yD^,^p{x - y)[Sl^{y) + .rj{y)\ (D12) 



for some propagator D^i,af} that will be determined shortly. Putting these solutions into (|D9j) gives the effective 
action 

Seff[{^K±}] = \jd'xj d^y [Sl\x) + r^{x)\ [Dl%,{x -y) + D^p.siv - ^)] [Sl\y) + J^'(y)] + OChl). 

(D13) 

Calculating the variation of the effective action using 

/i^^(i,x) = h+f,^(t,x) ^ ^ (D14) 



(D15) 



dVr(t,x) 

gives the gravitational waves radiated by the binary 

V(t,x) = 1 [ d^y[D:,%,ix-y)+D^f,^s{y-x)]S'''iy;^K). (D16) 



Self-consistency demands ()D16P to equal (jPlip when x^i = xx2 = and J!^'^ = 0. Therefore, Dap-ys{x — y) must 
be the advanced propagator. With this choice it is straightforward to show that (|D16[) is equivalent to (|4.28p . in the 
TT gauge, when the y-integral is evaluated. 

The fact that h- satisfies advanced boundary conditions is not a problem because x/f „ and Jj!" are set to zero at 
the end of the calculations (implying h- vanishes, from (jD12p ) so that there is no contribution to the equations of 
motion and gravitational waves from advanced radiation. The effective action is therefore 

5e//[{xA'±}] = Jd^x Jd''y[S'l^ix) + r^ix)]D:;^^,ix~y)[Sl\y) + Jl'{y)]+OChl). (D17) 
We can write ()D17P in terms of i times the matrix of propagators introduced earlier in the in-in formalism (j4.19p 

/ A T~)adv \ 

D'Jp,s = Jet ) (D18) 
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with a, & = ± so that (|D17p becomes 

Seff[{^K±}] = \j d^x J d'y[Sfix) + jf{x)]D^Xsi^-y)[S^,\y) + r,\y)] +0(P±). (D19) 

The CTP indices a, b are raised and lowered with the metric Cab in (|3.32[) . It is straightforward to show that setting 
the auxihary sources to zero in (jD19[) and varying the effective action with respect to x^- gives the radiation reaction 
force nil, m, derived above in (g^l). Indeed, ((DT91) is simply the in- in effective action of NRGR ([T^ with 
factors of i rearranged accordingly. 

The formulation of the radiation sector discussed here can be straightforwardly applied to the extreme mass ratio 
inspiral of a small compact object. In this context, there are only radiation gravitons and no potential gravitons to 
integrate out. Using an appropriate regularization scheme in the presence of a curved background spacetime we have 
calculated (bu t do not present the details here) the well-known self-force on the small compact object that was first 
derived in Indeed, the calculation is nearly identical (up to factors of i from applying the Feynman rules) to 

that given in [37|, which uses the in- in path integral framework. 
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